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Abstract 

It is known that if neutrino is a Dirac fermion with magnetic moment, then 
vl ^ ^ transition of supernova neutrinos may occur. The first stage of 
such transition is due to the neutrino spin flip inside the hot dense supernova core, 
while the second one - due to the neutrino spin precession in the galactic magnetic 
field on the way from the supernova to terrestrial detectors. This can result in 
the detection of 60-200 MeV neutrinos simultaneously with the "normal" supernova 
neutrino signal, which would be a smoking gun for the Dirac neutrino magnetic 
moment, [i. We argue that in case of a nearby supernova explosion (~ 10 kpc away 
from the Earth) one may observe such high-energy events in Super-Kamiokande if 
/X > 10~^^/iB, and in a Mt-scale detector if ^ > 0.5 • 10"^^//^- Such an observation 
by itself, however, may be not sufficient to determine the value of the magnetic 
moment, because of the ignorance of the interstellar magnetic field. We point out 
that if in addition a deficit of the neutronization burst neutrinos is established, it 
would be possible to extract the value of the magnetic moment from observations. 
We also briefly discuss a possible Majorana magnetic moment signature due to 
— > z^^t,r flip inside the supernova core. 
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1 Introduction 

It is natural for neutrinos to have magnetic moments. In the Standard Model (with 
massive Dirac neutrinos included) magnetic moment of the neutrino mass eigenstate I'i is 
proportional to its mass rrij and reads (see P,[2]) 
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In the extensions of the Standard Model neutrino magnetic moments may be orders 
of magnitude larger (see [S] and references therein). The best direct experimental limit 
currently reads |1] 

/i^, < 6 ■ 10-"/iB, 90%CL. (2) 

Astrophysical limits on the neutrino magnetic moment are two orders of magnitude 
better, though model dependent and less definite. Most of such limits come from con- 
sidering energy losses due to emission of neutrinos for various types of stars, such as red 
giants and white dwarfs, as well as for collapsing supernovae. Such energy losses, being 
dependent on magnetic moment /x, should not be too high. This allows to restrict neutrino 
magnetic moments (see [3] for the review and [5] for probably the most stringent limit of 
this type): 

/i < (1 - 3) ■ IO-^Vb. (3) 

Explosion of a nearby supernova, SN1987A, provided a possibility to constrain /i using 
non-observation of high-energy neutrinos {E^, = O(lOOMeV)) in two water Cherenkov 
detectors, which were operational at that time, 1MB and Kamiokande. Different aspects 
of this method were discussed and investigated in a number of papers [2], [6]-[ll]. The 
idea of the method is the following. Left electron neutrinos, created in the hot collapsing 
supernova core at the first second of the collapse, experience spin flips in collisions with 
charged particles, among which electrons and protons play the major role: 



vl+P ^ VR+p. (4) 

Sterile right neutrinos escape from the core without energy loss. Typical energies of this 
leaking neutrinos, O(lOOMeV), are of order of typical energies of the electrons inside the 
core (compare with 10-40 MeV for "standard" supernova neutrinos). Then neutrinos 
travel through the galaxy to reach detectors at the Earth, and a backward transition, 
— ^ t/^, may occur in the galactic magnetic field. If the neutrino magnetic moment is 
large enough, this would result in the detection of high-energy neutrinodj simultaneously 
with the "normal" supernova neutrino signal. As no such high-energy neutrinos were 
detected by Kamiokande and 1MB, one can constrain neutrino magnetic moment [TO] [11]: 

< IO-^Vb. (5) 

Here it is supposed that one half of right neutrinos turned to left neutrinos in the galactic 
magnetic field. 

In this letter we discuss prospects of exploring neutrino magnetic moment in the future 
collapse of a nearby Supernova. We concentrate mainly on the case of Dirac neutrinos, 
and only in section 3 briefly comment upon one possible Majorana neutrino signature. 
We estimate the expected number of high-energy neutrino events in Super-Kamiokande 

^ In this letter "high energy neutrinos" mean 60-200 Mev neutrinos. 



2 



(22.5 kt fiducial mass) and in a future Megaton-scale detector. It is found that, for a 
conservative estimate, Super-Kamiokande may be sensitive to /i > 10~^^fiB, while a 0.5 
Mt water detector - to fx > 0.5 ■ lO'^^/isO It is also pointed out that the deficit of 
neutronization burst neutrinos due to ul —>■ vr transition in the interstellar magnetic 
field may be another clear signature of Dirac neutrino magnetic moment. If both deficit 
of neutronization burst neutrinos and high energy neutrinos are observed, this may allow 
to determine the value of the neutrino magnetic moment. 

Two remarks should be made. First, in principle it is possible that the backward 
transition, ^l, occurs already in the star envelope, if a strong magnetic field is 

present there. In this case high-energy refiipped left neutrinos interact with the envelope 
matter and release their energy. This helps to produce a successful supernova explosion, 
but cancels the high energy neutrino burst. It was Dar [B] who first pointed out on this 
possibility. The effective matter potential, which is felt by left neutrinos and is not felt 
by right neutrinos, in general suppresses — > i^l transition in the magnetic field inside 
the star; nevertheless it may vanish at some radius, thus allowing z/r — > ul transition. 
This question was discussed in a number of papers [12]-[l5]. We suppose that ur ul 
transition does not occur inside the star. According to [T2], a sufficient condition for this 
is 



B being the magnetic field at ~ 100 km from the center of the star. 

The second remark concerns the exact meaning of /i. In general, electromagnetic prop- 
erties of neutrinos are described by 3 x 3 complex matrix M^m, which embraces diagonal 
(for Dirac neutrinos) and transition (both for Dirac and Majorana neutrinos) magnetic 
and electric dipole moments [16]. For ultra-relativistic neutrinos the phenomenological 
manifestations of electric dipole moments coincide with those of magnetic dipole moments 
|17j . Bearing this in mind, we do not consider explicitly electric dipole moments. Deal- 
ing with Dirac neutrinos, we assume for simplicity that M is diagonal in any basis, i.e. 
Mem = ■ 1, where 1 is a unit matrix and /i is a universal magnetic moment. 

The paper is organized as follows. In section 2 we obtain our main results. We pro- 
ceed in two steps: first we separately discuss neutrino propagation in the interstellar space 
(section 2.1), and then calculate the expected number of high-energy events in the de- 
tector (section 2.2). In section 3 we discuss the uncertainties involved and prospects for 
improving the accuracy of the results; also we make some remarks concerning diffuse neu- 
trinos from all past supernovae, as well as possible Majorana neutrino magnetic moment 
signature. The results are summarized in section 4. 

^ For n < 10~^^^B the number of events in the detector is proportional to /i"* (see section 2.1). This 
explains why increasing the fiducial mass by factor 20 results only in doubling of the sensitivity. 




(6) 



3 



2 Dirac neutrino magnetic moment signatures in the 
supernova signal 



2.1 Neutrino propagation in the interstellar space and vj^ 
transition in the galactic magnetic field 

The interaction of the neutrino magnetic moment with the magnetic field B leads to the 
neutrino spin precession (or, in other words, to ^ vl oscillations), which is described 
(in the ultra- relativistic case) by^l6j 

i^p{x) = {E + ij,(TB^{x))iy{x). (7) 



Here z/(x) = 




, E is neutrino energy, cr is a vector constructed from Pauli 



matrices, and Bx(x) is a component of B normal to the neutrino momentum. If for every 
X magnetic field Bx(x) lies in the same plane, then the phase of oscillations is given by 

= j ixB±{x)dx. (8) 

The oscillation probability, P{i'r ul) = sin^ 0, may be easily calculated for this case. 
For the constant magnetic field one gets 

,.,„B,..0,9(^)(^)(^), (9) 

Galactic magnetic field has a complicated structure (see, for example, [H], [19]). Its 
typical strength is not less than 1 /xG, and probably somewhat larger. It can be represented 
as the sum of regular (large-scale) and random (small-scale) components. Length scales 
of the random component are much smaller than 1 kpc, therefore, according to ([H]), this 
component is irrelevant for our purposes]! Length scales of the regular component are of 
order of 1 kpc. In the galactic disk regular magnetic field is directed along the spiral arms, 
clockwise or counterclockwise depending on the spiral arm. There is a variety of galactic 
magnetic field models (see the above mentioned references). We use a model described in 
[IH]. It reproduces the main qualitative features of the radial dependence of the magnetic 
field in the inner galaxy, i.e. two field reversals at ~ 4.5 and ~ 6.5 kpc from the center of 
the galaxy, as well as the characteristic strength of the magnetic field. 

We consider a frequently discussed case of a supernova exploding in the inner part of 
the disk of our galaxy, D = 10 kpc away from the Solar system. For simplicity, we assume 
that it is situated on the line which connects Solar system and galactic center. The radial 

^ Strictly speaking, relevance of the random magnetic field to the spin rotation is determined by 
7 — ji^ {B^)LcX [20], where Lc is a field length scale. Taking ji = 10~^^/xb, a; = 10 kpc and B ^ 1 /j,Gs, 
Lc ^ 10 pc (see [TH] and reference therein), one obtais 7 ~ 0.1 <C 1. This means that the effect of the 
random magnetic field on the spin precession may be neglected. 
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dependence of Bj_{r) reads [18] 



0.9 fiG kpc < r < 2 kpc 

3.8 fiG 2 kpc < r < 3 kpc 
3.1 /iG 3 kpc < r < 4 kpc 

BAr)= I -2.2 iiG 4 kpc < r < 5 kpc (10) 

— 1.9 jiG 5 kpc < r < 6 kpc 

1.9 /iG 6 kpc < r < 7 kpc 
2.5 yuG 7 kpc < r < 8 kpc, 

Here r is the galactocentric distance. The distance from the Sun to the galactic center is 
taken to be 7.2 kpc in [TH] . 

From (IHj) and ( fTOl) one obtains the probabihty of — * oscillations: 



sm 



1.1- 



/i 



(11) 



For /i > 10~^^/iB the probability oscillates rapidly with /i, phase being strongly de- 
pendent on the galactic magnetic field model. This means, in fact, that for such values of 
/i one should consider the phase as a uniformly distributed random value. In this case 



is also a random value. Its expectation value is Pav = 0.5, and Pi 



l^R-^l^L 



_ > 0.025 
with 90% probability. 

One should take into account flavour transformations along with the spin precession 
in the problem involved. Right electron neutrinos, produced in a supernova, u^r, quickly 
decohere into the the mixture of and i'2r (see, for example, [21] for numerical estimates; 
some subtle details of this phenomenon were investigated in [22] [23]). The corresponding 
fractions in the mixture are equal to cos^ 9i2 and sin^ Ou, Ou — 30° being the mixing angle. 
In the interstellar medium uir — > un and ^'2L transitions occur with the probability 

PuR^UL &s described above. Finally, in the detector vn and 1^21 show themselves as 
electron neutrinos with probabilities cos^ 6*12 and sin^ 612 correspondingly. Combining all 
the probabilities together, one finds the fraction k of right electron neutrinos, VeR-, which 
are converted to left electron neutrinos, VeL-, in the interstellar space on the way from the 
supernova to the Earth: 



(cos'^ 9i2 + sin^ 9i2, 



P 



(1 -0.5sin2 2^i2)P^ 



0.6P, 



^R^VL ■ 



(12) 



2.2 Dirac neutrino magnetic moment signatures in the detectors 

In ref. [IDj an expression for the number of events with energies greater than 60 MeV in 
the water Cherenkov detector is suggestecfl: 

---'° °«-«)-^--- -^ (10^)^ (^) (lo^)"^- 

^To be precise, our ea. p3p may be obtained from eq.(ll) in ref. lOj, if one introduces conversion 
probability n in their expression and conservatively chooses the emission time to be 0.5 s (it is 
considered to be (0.5-l)s in lOJ) 
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fiducial mass 
of H2O 


number of events, 


number of events, 

U < W~^^Un 


Super-Kamiokande 


22.5 kt 


1.4 PuR^UL ( 10-13 /xs^ 


I 


<1 


MEMPHYS 


440 kt 


26 PuR-*UL (10-13^5) 


'2 


22 ( 10-13 



Table 1: Estimated numbers of events in Super-Kamiokande and MEMPHYS for the 
supernova in the inner part of the galaxy, 10 kpc away from the Solar system. In case 
/i > 10~^^ fiB the present knowledge of the galactic magnetic field does not allow to obtain 
any reliable value for the phase of spin precession 0; a priory the probability of — >■ ul 
conversion, P^^^^^^ = sin^0, may take any value between and 1. 

The two-orders uncertainty mainly comes from the approximate treatment of state of 
matter in the hot supernova core. We postpone the discussion of this point to the next 
section, and conservatively use the lower value of N. The estimated numbers of events for 
Super-Kamiokande and for the proposed ~ 0.5 Mt detector MEMPHYS [240 are presented 
in Table 1. To obtain this estimate in case fi < 10~^^ fiB we substitute sine by its argument 
in eq. lfTT]) . 

As was mentioned in section 2.1, it is extremely hard to estimate the probability of 
transition, P^^^^j^, in case when yu > IQ^^^jiB- However, it is interesting to note 
that a supernova explosion may provide an independent possibility to measure P„^^„^. 
Namely, it was argued in [27j that the number of neutrinos from the neutronization burst 
is almost a model-independent value; the uncertainty in the neutrino event number at a 
Mt-scale detector was claimed to be (10-15)% for a supernova at 10 kpc. Therefore, if a 
considerable deficit of the neutronization burst neutrinos is revealed, it may be attributed 
to the Vl Ufi transition of this neutrinos due to their precession in the interstellar 
magnetic field. This by itself may be considered as an evidence for the Dirac magnetic 
moment of neutrino, though not a completely convincing one. The deficit is determined 
by i^L transition probabihty Pyj^^^^, which is equal to Py^^y,^. Thus one can extract 

this probability from the value of the deficit. If the high-energy neutrinos from the same 
explosion are detected, not only the evidence for the Dirac magnetic moment of neutrino 
would be compelling, but it would be possible to determine fi from this observation, using 
a more precise version of eq.( fT3i) . 

However, it is worth noting that in some cases high-energy neutrino events from the 
supernova explosion may be observed, but the deficit of the neutronization burst neutrinos 
may not be clearly established, and visa versa. For example, the first situation may take 
place if at the moment of a supernova explosion no Mt-scale detector is operational. In 
this case only a few neutronization burst neutrinos can be detected by Super-Kamiokande 
and other detectors of the similar sensitivity. An opposite possibility may be realised if a 
Mt-scale detector registers neutrinos from a supernova which possesses a strong magnetic 
field. In this case condition ([U]) may be violated, and one may observe the deficit of the 

^There are other Mt-scale water detector proposals, e.g. Hvper-Kamiokande ^5] and UNOpS]. 
Planned sensitivities of all such detectors are comparable. 
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neutronization burst neutrinos without detecting high energy neutrinos. 



3 Discussion 

First of all, let us discuss eq.( |T3|) . It is obtained from 
N = K- I dE 



Mh^o (^{E) f ,^ f ,3„dn^R 



mHaO 4:7!'D'^ J J dEdt 



/ rf'r-^(E,p(r,t),T(r,t)) 



(14) 



Here 

K is the fraction of v^r converted to v^l in the interstellar space on the way from the 
supernova to the Earth (see section 2.1), 

MhjO and mnjO are the fiducial mass of water and the mass of the H2O molecule, corre- 
spondingly, 

D is the distance from the supernova, 

(j(£') is the cross section of the reaction z/g + O e~ + F for the neutrino with energy E, 
dn^^{E, p(r,t),T{r,t))/dEdt is the number of rigt neutrinos with energy E emitted per 
unit energy interval in unite time from unite volume of supernova matter with density 
p{r,t) and temperature T{r,t). This later quantity is integrated over the volume of the 
the supernova core and over the time of explosion. Two nontrivial points are to calculate 
dn^^/dEdt for given E,T,p, and to perform the integration over d^rdt. 

The first point mainly deals with the calculation of probabilities of the neutrino spin 
flips on electrons and protons (reactions (jlj)) in a dense hot media. In [TU] the media 
effects were accounted for by introducing the effective mass in the photon propagator. 
More sophisticated considerations [5] [28] give a somewhat larger result for driy^/ dEdt. It 
is not completely clear how the spin flip on nuclei contributes during different stages of 
supernova explosion (for example, Barbieri and MohapatrapU] consider only spin flips on 
protons (and electrons) after the bounce, while Notzold[llJ relies completely on the spin 
flips on large nuclei before the bounce). 

To proceed with the second point, one, in general, needs to know the evolution of the 
state of matter during the explosion, i.e. p{r,t) and T{r,t). In ref. [10] an oversimplified 
model was used: both temperature and density were considered to be constant within a 
sphere of volume V (supernova core) for a time interval t. The following numerical values 
were used: 

V = A ■ 10^^ cm^ (which corresponds to the core radius 10 km), 

t=(0.5-l) s, 

p = 8 ■ 10^^ g/cm^, 

T = (30 - 60) MeV. 

From the modern point of view (see, for example, [22]) typical density and temperature 
are overestimated in [10]. This could lead to the overestimation of the of the z/^j production 
rate. We are not aware of any calculations of this quantity which accurately account for 
the evolution of the state of the supernova core. Such calculations could greatly diminish 
the uncertainties in the predicted number of high energy neutrino events. 

We think that underestimation of the spin flip probability and overestimation of the 
supernova core temperature and density in [TU| could cancel each other to some extent. 
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Nevertheless, in view of this uncertainties we used the lower value for from the interval 
suggested in [ TD1 |. 

We would like to point out that the detailed calculations of neutrino spin flip rate 
with supernova dynamics included are of interest not only in the case of Dirac neutrinos, 
but in the case of Majorana neutrinos also. Majorana neutrinos may have transition 
magnetic moments, which can lead to Ve ^fi,T spin flips incide the core. Non-electron 
antineutrinos, and Vri may be then partially converted to electron antineutrinos, i>e, 
due to the flavor transformations. As in the first tens of milliseconds from the beginning 
of the collapse only electron neutrinos are expected to be produceci§ (see, for example, 
[27] [29]), detection of electron antineutrinos at this stage would be a smoking gun for the 
Majorana neutrino magnetic momentJil From estimates presented in ref. [5] one can get 
that for /i = 10~^^/iB and for the oversimplified supernova core model, described above, 
with T = 30 MeV, energy output of ~ 150 MeV flipped neutrinos is of order of lO^^erg/30 
ms. In the case of Majorana neutrinos this gives ~ 10^°erg/30 ms in cooled (~ 10 MeV) 
non-electron antineutrinos. This is only one order of magnitude less than the energy 
output of the neutronization burst (which duration is (20-30) ms). However one should 
carefully take into account supernova dynamics, as well as antineutrino diffusion from the 
central part of the core to the neutrino sphere, in order to justify (or refute) this naive 
estimate. 

There is an ongoing activity for detecting the Diffuse Supernova Neutrino Background 
(DSNB), i.e. neutrinos emitted by all previously exploded supernovae in the universe 
(see [33] for the most stringent experimental limit and ^34j for the theoretical review). 
A question may arise whether the high-energy supernova neutrinos form a detectable 
analogue of the DSNB. Unfortunately, simple estimates show that this is unlikely. A 
rough estimate for the event rate in the detector, dN/dt, may be obtained: 

^ ~ 47r(10kpc)2/?oiV^, (15) 
dt H 

where c is the speed of light, H = 70 km/s Mpc~^ is the Hubble constant, Rq ~ 10"'* yr~^ 
Mpc~^ is the local supernova rate, and N is defined above (see eq. ffT^ ). For = lO^^^/x^ 
one gets dN/dt ~ 1 yr"^ in the MEMPHYS detector. It is much smaller than the expected 
rate of atmospheric neutrino events with energies less than 200 MeV, which is of order of 
1000 yr-i[35]. 

4 Summary 

Neutrino signal from a future nearby supernova is likely to allow to probe Dirac neutrino 
magnetic moment /i up to /i > 10~^^/is at Super-Kamiokande and /i > 0.5 • lO^^^/i^ 
at a future Mt-scale detector. A smoking gun for such values of /i is the detection of 
(60-200) MeV neutrinos. With the Mt-scale detector another clear signature of the Dirac 

^In particular, this is the case for the beginning of the neutronization burst. 

^This effect does not depend on the supernova magnetic field strength, and therefore differs from the 
similar one based on the I'e —^ v^^t transition in the presence of the strong supernova magnetic field (see 
[30]- [H2] and references therein). 
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magnetic moment may be observed: the deficit of tlie neutronization burst neutrinos. If 
botli signatures are observed, one is able to extract the value of the magnetic moment 
from the observations. 
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